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ABSTRACT 


In outer space , the only means of rejecting heat is that of thermal radia- 
tion . For this purpose heat rejection space radiators have been developed. 
Internally generated heat within the space craft is conveyed to the radiator 
by a circulating coolant (e.g.,a water-glycol mixture). The coolant circulates 
in tubes (usually in a parallel circuit ) which are joined by finned surfaces. 
Heat from the coolant conducts into the fin and is dissipated by the radiation 
to the surroundings . The finned surface may be subjected to an external 
irradiation firom a source such as the sun . 

In this work , a finite- difference solution is obtained for ID radiating 
fins of rectangular, parabolic and trapezoidal profile, the material being C- 
steel, aluminium and copper. It is also assumed that the conductivity is 
a function of temperature . The results for the limiting case of constant 
conductivity rectangular fin have been checked with the earlier solution of 
NASA and found very satisfactory. Based on simulation results, that is , the 
fin temperature distribution and the heat transfer , an optimization study 
has been carried out with regard to the fin geometry and dimensions as well 

vii 



as materials in order to maximize heat transfer from the fin. An exhaustive 
search method has been used. It is found that aluminium is the most suitable 
material. Also, trapezoidal and parabolic fins are recommended for efficient 
heat transfer. 
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Nomenclature 


a 

b 

c 


Cl 


C2 


G 

k 

ki 

T 

Ti 


profile of parabolic fin is y = a + 

profile of parabolic fin is y = a + + cx^ 

profile of parabolic fin is y = a + + cx^ 

— _EL_ 

_ 1-0 
creLH^ 

external irradiation 

thermal conductivity of material 

is the thermal conductivity of material at grid point 

non dimensional temperature 

non dimensional temperature at i*'* grid point 

non dimensional temperature at grid point in previous itera- 
tion 
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to base temperature of the fin 

ti =— 

s ae 

V volume of the fin 

W width of rectangular fin 

Wi width of trapezoidal fin at base 

W 2 width of trapezoidal fin at tip 

X non dimensional distance toward the fin 

X distance 

Greek letters 

a absorptivity 

e emissivity 

a Stefan - Boltzmann’s constant 5.67 x 10~^W/m^.K^ 

p density of material 


X 



Chapter 1 


INTRODUCTION 

1.1 GENERAL REMARKS 

One of the most important problems of modern space technology is that 
of thermal control-the maintenance of a desired temperature level with in 
space craft (manned or unmanned) or satellite . Such spacecraft subjected 
to heat input of both internal and external sources . Internal sources consist 
of electrically generated heat , while the external sources consist of thermal 
irradiation from the sun , planetary bodies ,etc . If a net gain of heat is 
experienced, some means of heat rejection must be provided in order that 
thermal control be maintained. 

In outer space , the only means of rejecting heat is that of thermal radi- 
ation . For this purpose heat rejection space radiators have been developed. 
A typical space radiator configuration is shown in Fig 1.1. 
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Figure 1.1; Heat rejection space radiator 

Internally generated heat with in space craft is conveyed to the radiator 
by a circulating coolant (e.g., a water glycol mixture ) . The coolant circu- 
lates in tubes (usually connected in parallel circuit) which are joined by the 
finned surfaces . Heat from the coolant conducts into the fin and dissipated 
by radiation to the surroundings. The finned surface may subjected to an 
external irradiation from a source such as sun .The basic question to be an- 
swered is: Given a radiating fin of known geometry, maintained at a base 
(i.e.tube temperature) and subjected to external radiation, what is the rate 
at which heat is dissipated by the fin to the environment? 

1.2 LITERATURE REVIEW 

Optimum design of radiating fins has been the subject of research by a num- 
ber of investigators. Wilkins’ (1960) analysis of a single fin freely radiating 
to space shows that the profile of the fin of least material follows a power 
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lawji.e.y , where y is the fin thickness at the length x from the tip of the 
fin and a = 3.5 . Sparrow et al. (1961) presented the optimum design of rar 
diatively interacting longitudinal fins without considering fin -to -base mutual 
irradiation. Later Sparrow and Eckert (1962) emphasized the importance of 
mutual radiation interaction between the fin and its base surface. The work 
by Kerlekar and Chao (1963) also neglected fin-to-base interaction while ana- 
lyzing the optimum design of trapezoidal fin. Schnurr et al. (1976) employed 
a nonlinear optimization technique to determine the minimum weight design 
for straight and circular fins of rectangular and triangular profiles, protruding 
from a cylinder, considering fin-to-fin and fin-to-base radiation interactions 
. Chung and Zhang (1991a) determined the optimum shape and minimum 
mass of a thin fin accounting for fin - to - base interaction based on a vari- 
ational calculus approach . Chung and Zhang (1991b) later extended their 
analysis to minimize the weight of a radiating fin array projecting from a 
cylindrical surface considering both fin-to-fin and fin to base interactions . 
The results are presented in form of graphs and no functional form is given for 
the optimum profile of the fin. Krishnaprakas (1996) presented the optimum 
design of a rectangular plate fin array extending from a plane wall employing 
a nonlinear optimization method . Krishnaprakas (1997) presented the opti- 
mum length of a radiating longitudinal fin array extending from a cylindrical 
surface. Chung et al. (1996) analyzed and optimized a new designed for a 
four-fin radiating array in a fuzzy environment. Here, the objective function 
is not precisely defined , rather the two conflicting goals of weight reduction 
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and horizontal height minimization for space saving are fuzzily defined as 
fairly light and reasonably short. 

1.3 OBJECTIVES 

Radiating fins play an important role in thermal control of space vehicles and 
satellites. Since the weight and material are the major design considerations 
in heat exchangers , it is highly desirable to obtain an optimum design with 
respect to fin mass. The objectives of this study are to numerically simulate 
heat transfer from radiating fins of rectangular, parabolic and trapezoidal 
profiles and find the dimensions of the fins of different materials for maximum 
heat transfer. 
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Chapter 2 


PROBLEM FORMULATION 
AND METHOD OF SOLUTION 

2.1 PROBLEM STATEMENT 

If the fluid in adjacent radiator tubes is at same temperature, the fin may- 
represented by the models in the Fig2.1,Fig2.2 and Fig2.3. The base tem- 
perature of the fin is to . The other end of the fin has the condition ^ = 0.0 
. The fin material has the thermal conductivity k which may be function 
of temperature, and expose surface has total emissivity e . The surface is 
exposed to an external irradiation G (perhaps solar) to which it exhibits a 
total absorptivity a, since such radiators are often made an integral part of 
space craft skin, only one of the surface of fin is considered to be radiatively 
active . 
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2.2 GOVERNING DIFFERENTIAL EQUATIONS 


AND BOUNDARY CONDITIONS 


2.2.1 FOR RECTANGULAR FIN 

If the conduction in the is fin taken to be one-dimensional, then the heat 
balance taken on an element 6x in length yields the following relation (for 
a unit depth) of the rectangular fin of length L and uniform thickness W 
(Fig.2.1) 


-kW^ |a:= -kW^ \x+5x +ecrt^Sx - aG5x (2.1) 

CLOO CLX 

’A: ’ may be also function of temperature. The actual functions of various 
fin materials used in this study are shown in Appendix ’B’. 

No convective loss at the fin surface is included. For 5x 0.0 , the 
following differential equation must be satisfied for each point. 


d /tfirdt-. /,4 cxG .,^ 

—(kW—) = ear{t ) 

dx^ dx^ €cr 


( 2 . 2 ) 


The non-homogeneous term is recognized as the equilibrium temperature, 
or equivalent sink temperature, defined in equation. This temperature is the 
equilibrium temperature the surface would achieve if isolated in irradiation 
G. Thus with 
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Figure 2.1: RECTANGULAR FIN 
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(2,3) 


oG 


ea 


The differential equation for the temperature distribution in the fin is 


- ti) = 0.0 (2.4) 

Boundary conditions are 
At x = 0.0 


At X = L 


t = to 


(2.5) 


dt 

dx 


0.0 


( 2 . 6 ) 


The solution of above equation is best discussed in terms of the following 
dimensionless variables. 


X = 


X 

L 


T = 


L 

to 



8 



Introduction of these definitions in eq. (2.4) yields 


^)-0.0 (2.7) 

The dimensionless boundary conditions to be satisfied are 
AtX = 0.0 


T = 1.0 


( 2 . 8 ) 


At X = 1.0 


^ = 0.0 (2.9) 

2.2.2 FOR TRAPEZOIDAL FIN 

If the conduction in the fin is taken to be one-dimensional, then the heat 
balance taken on an element Sx in length yields the following relation (for a 
unit depth) for the trapezoidal fin of length L .The fin width at a: = 0.0 is 
Wimdx = LisW 2 . See Fig. 2.2. 

’k ’ may be also function of temperature. The actual functions of various 
fin materials used in this study are shown in Appendix ’B’. 


-k{Wi-iWi-Wi)j)^ |*= -k(Wi-iWi-W 2 )^)^ U+Sx +€atHx-aGSx 

( 2 . 10 ) 
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No convective loss at the fin surface is included. For 5x 0.0 , the 
following differential equation must be satisfied for each point. 

The non-homogeneous term is recognized as the equilibrium temperature, 
or equivalent sink temperature, defined in equation. This temperature is the 
equilibrium temperature the surface would achieve if isolation in irradiation 
G. Thus with 


aG 


ea 


( 2 . 12 ) 


The differential equation for the temperature distribution in the fin is 


- (Wi - W2)~)^) - - tt) = 0.0 (2.13) 

ax L dx 

Boundary conditions are 
At x = 0.0 


t = to 


(2.14) 


At X — L 


^ = 0.0 (2.15) 

Qi2C 

The solution of above equation equation is best discussed in terms of the 
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following dimensionless variables. 



T = 


£ 

to 



Introduction of these definitions in eq.(2.13) yields 


- (r* - I?) = 0.0 (2.16) 

The dimensionless boundary conditions to be satisfied are 
At X = 0.0 


T = 1.0 


At i: = 1.0 



(2.17) 


(2.18) 
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2.2.3 FOR PARABOLIC FIN 


If the conduction in the fin is taken to be one-dimensional, then the heat 
balance taken on an element 5x in length yields the following relation (for a 
unit depth) for the parabolic fin of length L .If origin is taken at a; = 0.0 and 
at the middle point of width and positive x axis is taken in the direction of 
fin length . The fin profile is of second degree parabola . See Fig.2.3. 

'k ’ may be also function of temperature. The actual functions of various 
fin materials used in this study are shown in Appendix ’B’. 

Upper part of fin parabola 

y = {a + bx + crc^) 

Lower part of fin parabola 

y = -(a + bx + cx^) 


~-2k{a+hx+cx'^)^ la:= -2k{a-{-bx-\-cx‘^)~ \x+6x +eat'^dx-aG6x (2.19) 

No convective loss at the fin surface is included. For 5x 0.0 , the 
following differential equation must be satisfied for each point. 
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( 2 . 20 ) 


-^{2k{a + bx + cx^)^) = ecr(i^ — 
dx' " 'dx^ ^ ea^ 

The non-homogeneous term is recognized as the equilibrium temperature, 
or equivalent sink temperature, defined in equation. This temperature is the 
equilibrium temperature the surface would achieve if isolation in irradiation 
G. Thus with 


e<T 


( 2 . 21 ) 


The differential equation for the temperature distribution in the fin is 


■^(2k(a + bx + cx^)^) - €cr{t^ — t,) = 0.0 (2.22) 

dx dx 

Boundary conditions are 

At X = 0.0 


t = to 


(2.23) 


At X = L 


^ = 0.0 (2.24) 

ax 

The solution of above equation equation is best discussed in terms of the 
following dimensionless variable. 
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r = 


£ 

to 


T,= 


h. 

to 


Introduction of these definitions in eq. (2.22) yields 


- (T- - 1?) = 0.0 (2.25) 

The boundary conditions to be satisfied are 
At ;«C = 0.0 

T = 1.0 (2.26) 

At A = 1.0 

% = 0.0 (2.27) 


2.3 METHOD OF SOLUTION 


The finite-difference method is used to discretize the governing differential 
equations. The finite difference method enables one to integrate a differential 
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equation numerically by evaluating the values of the function at a discrete 
(finite) number of points . The origin of this method is Taylor series ex- 
pansion which assumes that the function is smooth , that is continuous and 
differentiable, using central difference 


dT 

dX 


2AX 


+ 0{AXf 


Ti^i 2Ti-\- Ti—i r\( w\^ 

(Axp — + 


(2.28) 

(2.29) 


The notation 0{AXY means that in arriving above, term of order of 
(AX)^ and higher have been neglected. 0(AX)^ is called the truncation 
error. The truncation error is the difference between the exact mathematical 


expression and numerical approximation. 


2.3.1 DISCRETIZATION FOR RECTANGULAR FIN 


Assuming 



(2.30) 


Governing equation 


= (2.31) 

Boundary conditions 
At A = 0.0 . 
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Figure 2.4: Grid point for Rectangular fin 


T = 1.0 (2.32) 

AtX = L 

§=0.0 (2.33) 

The above equation can be written as 

To indicate non-dinaensional temperature in previous iteration. 

For i = 1 


18 



Ti = 1.0 


(2.35) 


ki is the thermal conductivity of material at point at Toi . 
For i = 2 


+ 1.0 - 2Ti Ti+i 

Oi/Ci h Cl( 


AX2 


2AX 




ki-i 


2AX 


) - (T^-T, - r/) = 0.0 
(2.36) 


Tii- 


AX^ 


■^oiJ 




= -C'x( 


4ki 


Ci+l 


+ k 


4AA'2 


-)-T’,^ 


(2.37) 


For i = 3, , n — 2, n — 1 


Ciki 


Ti+i + Ti-i — 2^ ^ ^ — Ti_ 


AX2 


2AX 


4( 


«t+i 


ki- 


2AX 


-) - (TSiT, - Jj) = 0.0 


(2.38) 




— ki+i + ki-i). - 2^71 A:,- ^3\ , , C'i(A:i4.i - ki-i), 

4AX2 AX2 2AX2 ^ 



(2.39) 


For i = n using image-point technique, eq.(2.39) takes the form 


19 



rp \ , rp / ‘ICiki rp3\ rp4 

Iv, A y 2 / A V2 ^ 0 %) "^5 


(2.40) 


2.3.2 DISCRETIZATION FOR TRAPEZOIDAL FIN 



Figure 2.5: Grid points for Trapezoidal fin 


Assuming 


_ 1.0 

ecrtlL^ 


(2.41) 




-T?Ti - Tf = 0.0 


(2.42) 
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Assuming 


Si = Wi- {Wi - W2)Xi (2.43) 

To indicates non-dimensional temperature in previous iteration. 

For i = 1 


Ti = 1.0 


(2.44) 


ki is the thermal conductivity of material at point at 
For i = 2 

- kiC2{Wr - » 2 ) 




2AX 

2C2kiSi ^3,^ ^ ,SiC2ki C2ki{Wi — W 2 ) ^ o. (j^i+i ~ 


(2.45) 


,SiC2ki , ^ W-W 2 ) 
^ AA2 ^ 2AA 


-(725i 


fct+i ~ fej-i ’' 
4AA2 


(2.46) 


For i = 3, n — 2, n — 1 


+CjS,( 


T, 


1 + 1 


Ti_i^ — ki-x 


2AA 


-)( 


2AA 


) - {T%Ti 


Tf) 


0.0 


(2.47) 
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^ CMWi -W2) , ^ ^ {ki+r - 

AX2 ^oiJ + ii+i(-^ 2AX +C^'5i . ^^^---) 


1 71 Ai^^ki „ fci(W^i W2) _ ^ fcj+i — /ij-i 4 

+Jz-ll-^j2- + ^2 — ;a V2 ) = -^5 


2AX 


4AX2 


(2.48) 


For i = n , using image-point technique, 


Ti-i = Ti 


t+i 


(2.49) 


Putting eq.(2.49) in eq. (2.48), 


_ .^ CAij^s . rp / "^ C^Siki ^ — 3 . _' T ’4 

■‘i— H ^ y 2 /"T'-^iV A V2 "^oi/ s 


AX2 


(2.50) 


2.3.3 DISCRETIZATION FOR PARABOLIC FIN 


Assuming 


C2 = 


1.0 


eaL^I^ 


(2.51) 


Governing equation 


C2^ik2{a + bXL + CX'^L'^)—) - (T^ - T/) = 0.0 


(2.52) 


Boundary conditions 
At X = 0.0 . 
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Figure 2.6: Grid points for Parabolic fin 


r = i.o 


(2.53) 


AtX = L 


^ = 0.0 (2.54) 

dX 

The above equation can be written as 

C2k^{a+bXL+CX'^L'^)+^C2{a-^-hXL+CX‘^L^)^HhL+2XL'^)C2k^ 

-{T^Ti - r/) = 0.0 (2.55) 


Assuming 
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Sli = 2{a + bXL + cX^L‘^) 


(2.56) 


To indicates non-dimensional temperature in previous iteration . 
For i = 1 


Ti = 1.0 


(2.57) 


ki is the thermal conductivity of material at point at Tgi . 

For i = 2 

+ kiC2{hL + 2cJ^L2)S±^ 

- Mr, - J?) = 0,0 (2.58) 

2C2kiSli ^3, , ^ ,SUC2ki , C2ki{hL + 2cXL2) ^ 

2AX +C7251i ) 


^4 rSUC2ki ^ ki{hX + 2cXL2) ^ ki^, - ki^i, 

^ 2AX 4AX2 ^ 


(2.59) 


For i = 3, ,n -- 2,n - 1 

+ kiCiibL + 2cXL") ^^^i 


3 <(- 


+C25i( 

2C72A:i51i 

AX2 


t t-|-X 


ti-1 


2AX 


)( 


A^i+i 


fci. 


2AX 


i)-( 23 T,-I?) = 0.0 


,,,3, , C3*,(5L + 2cXL^) 

•*oij+-^HU /\ y 2 ' 9/\ y hOaOlj 


(2.60) 


(^i+1 ~ fej-l) - 
4AX2 ' 
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I 'T' S\{C2^i ^ hi(bX “1“ ^cXL ) ^ ^ ^i—l\ rp4 (cf ’\ 

^Ti-A-^-C2 2^3^ C2S\,-^^^)--T, (2.61) 

For i = n , using inaage-point technique , eq.(2.61) takes form 

Ti.,(?^^)+r,(-?^^-Ta = -T,'‘ ( 2 - 62 ) 

2.4 CALCULATION OF HEAT TRANSFER 
AND EFFICIENCY OF THE FIN 

The complete solution of the nonlinear equation here can be obtained the by 
an iterative method. Rather than seek the distribution , the principal item 
of interest is the heat dissipated by the fin . 

9=-W(^)a==o.o (2-63) 

W in case of trapezoidal fin 


W = Wi 


(2.64) 


W in case of parabolic fin of profile 
For upper part 
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y = {a + hx-\- cx^) 


For lower part 


y = —(a + bx + cx^) 


W = 2a (2.65) 

In form of non dimensionless variables the q 


Q = 


hWto 

L 


& 


( 2 . 66 ) 


The heat dissipated is best expressed in term of fin efiiciency, which is the 
ratio of q to the heat that would be dissipated if the whole fin is maintained 
at base temperature. 


effidmcu = 

Eq.(2.67) can also be written in other form also 
effide^ = _ rfl 


(2.67) 


( 2 . 68 ) 
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Where 



_ ecrtlL^ 
~ kW 


(2.70) 


2.5 METHOD OF SOLUTION (along with flow 
chart) 

The r^ulting equation after discretization are solved by the well known Tri- 
diagonal matrix algorithm (TDMA). The flow chart of the overall solution 
procedure is shown in Fig.2.7. A grid independence test has been performed 
for each type of fin. It is found that 20 grid points are sufficient. 
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Chapter 3 


RESULTS AND DISCUSSION 


Here the result are given in mainly two parts, that is, for constant conduc- 
tivity and variable conductivity. For all the results presented here, e = 0.5, 
a = 0.1, G = 315.4-2^. 

ftl 

3. 1 CONSTANT CONDUCTIVITY OF FIN MA- 
TERIALS 

3.1.1 RECTANGULAR FIN 

TEMPERATURE DISTRIBUTION 

The temperature of rectangular fin of steel of constant conductivity (A: = 

45 W/mK).See Fig. 3.1 
Taking 
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W = 0.000375m 


L = 0.075m 

COMPARISON WITH ’NASA’ RESULT 

The Efficiency Vs \/A for rectangular is compared with the S.Lieblein,NASA 
Tech. Note D-196, National Aeronautics and Space Administration , Wash- 
ington, D.C.1959. and is found very satisfactory (Fig.3.3). It is to be noted 
that the NASA results are available only in the graphical form in Chapman 
(1989) . The NASA results axe shown in Fig.3.2. 

3.1.2 TRAPEZOIDAL FIN 

TEMPERATURE DISTRIBUTION 

Temperature of trapezoidal fin of steel of constant conductivity (k = 45 
W/mK). See Fig. 3.4. 

Using 


Wi = 0.000375m 


W2 = 


Wi 

2.0 
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Temperature distribution T VS X for rectangular profile 



ejnjBjediue) iBUOjSueuiipuoN 

Figure 3.1: Temperature distribution of rectangular fin of constant conduc- 
tivity 


31 


Nondimensional distance 







to 

CO 



Aou 0 ! O |||3 

Figure 3.3: Efficiency Vs y/X for rectangular fin of constant conductivity 
*0.0’ indicates ^ = 0.0, ’0,5 ’indicates ^ = 0.5, ’0.7’ indicates ^ = 0.7 
’0.9 ’indicates ^ = 0.9 , ’NASA’ indicates original NASA result 
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SQRT(lamcla) 



L = 0.1125m 


EFFICIENCY VS sfl 
Here IV2 = ^ 

Here A is 


_ eat^L^ 
~ kWi 


(3.1) 


See Fig. 3.5. 


3.1.3 PARABOLIC FIN 

TEMPERATURE DISTRIBUTION 

For a parabolic fin, the graph of T Vs X is given here for steel of A: = 45 
W/mK. See fig. 3.6. 

Using 


a = 0.0001875 


b = -0.001125 
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Figure 3.4: T Vs X for Trapezoidal fin for constant conductivity 
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Aouapi^a uy 

Figure 3.5: Efficiency Vs y/X for Trapezoidal fin of constant conductivity 
’0.0’ indicates ^ = 0.0, ’0.5 ’indicates ^ = 0.5, ’0.7’ indicates ^ = 0.7 
’0.9 ’indicates ^ = 0.9 , 
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SQRT(LAMDA) 



c = 0.000005 


L = 0.1m 

EFFICIENCY VS VX 
Using 

a = 0.0001875 

b = -0.001125 

c = 0.000005 


L = 0.1m 


Here A is 


ecrtlL^ 

2ka 


(3.2) 
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Figure 3.6: T Vs X for parabolic fin for constant conductivity 
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0.84 





See Fig. 3.7. 


3.2 VARIABLE CONDUCTIVITY OF FIN MA- 
TERIAL 

3.2.1 RECTANGULAR FIN 

TEMPERATURE DISTRIBUTION 

The temperature distribution of Rectangular fin is drawn taking copper 
as the fin material. If the base temperature is low, then there is almost no 
variation in temperature distribution with respect to constant and variable 
conductivity . Therefore ,the base temperature is taken as 673K. See Fig. 
3.8. 

Talcing 


W = 0.000375m 


L = 0.075m 


EFFICIENCY VS ^/X 

Base temperature is taken as 323K. 

For carbon steel 
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Figure 3.7: Efficiency Vs y/\ for parabolic fin of constant conductivity 
^0 indicates ^ = 0.0, ’0.5 ’indicates ^ = 0.5, ’0.7’ indicates = 0.7 

’0.9 ’indicates = 0.9 , 


lO 

o 
o 
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SQRT(lamda) 




ajniBjadiua) fBuosuaoijpuoisj 

Figure 3.8; Comparison of T Vs X for constant conductivity and variable 
conductivity for rectangular fin of copper 

’c* indicates constant conductivity, ’k’ indicates variable conductivity 
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Nondimensonal distance 



_ ecrtlL^ 
~ kiW 


(3.3) 


ki is Thermal conductivity at base temperature. 

See Fig. 3.9. 

3.2.2 TRAPEZOIDAL FIN 

TEMPERATURE DISTRIBUTION 

The temperature distribution of Trapezoidal fin is drawn talcing copper 
as the fin material. If the base temperature is low, then there is almost no 
variation in temperature distribution with respect to constant and variable 
conductivity . Therefore ,the base temperature is taken as 673K. See Fig. 
3.10 . 

Using 


Wi = 0.000375m 


W2 = 


Wi 

2.0 


L = 0.1125m 


EFFICIENCY VS VX 
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CDOOr^ COIO'^COCNJ^ 

ooooooooo 


AouejOj^g 

Figure 3.9: Comparison of Efficiency Vs y/X for constant conductivity and 
variable conductivity for rectangular fin of carbon steel 
’0.0’ indicates ^ = 0.0, ’0.5 ’indicates ^ = 0.5, ’0.7’ indicates ^ = 0.7 
’0.9 ’indicates = 0.9 , 

For variable conductivity material ’O.Ov’ indicates ^ = 0.0, 

’0.5v’ indicates ^ = 0.5, ’0.7v’ indicates = 0.7, ’0.9v’ indicates ^ = 0.9. 
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Figure 3.10: Comparison of T Vs X for constant and variable thermal con- 
ductivity for trapezoidal fin of copper 

’clkr indicates constant conductivity. ’c2k2’ indicates variable conductivity 
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For caxbon steel for a particular case = ^ 
Base temperature is taken as 323K 


~ hWi 


(3.4) 


ki is the the thermal conductivity of carbon steel at base temperature. 
See Fig. 3.11 . 


3.2.3 PARABOLIC FIN 

TEMPERATURE DISTRIBUTION 

The temperature distribution of Parabolic fin is drawn taking copper as 
the fin material. If the base temperature is low, then there is almost no 
variation in temperature distribution with respect to constant and variable 
conductivity . Therefore ,the base temperature is taken as 673K. See Fig. 
3.12 . 

Using 


a = 0.0001875 


b = -0.001125 


c = 0.000005 
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Aou0jo|||3 

Figure 3.11: Comparison of Efficiency Vs\/A for constant conductivity and 

variable conductivity of trapezoidal of copper 

’0.0’ indicates ^ = 0.0, ’0.5 ’indicates ^ = 0.5, ’0.7’ indicates ^ = 0.7 

’0.9 ’indicates = 0.9 , 

For variable conductivity material ’O.Ov’ indicates ^ = 0.0, 

’0.5v’ indicates ^ = 0.5, ’0.7v’ indicates ^ = 0.7, ’0.9v’ indicates ^ = 0.9 
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EFFICIENCY VS 
Using 


L = 0.1m 


a = 0.0001875 


6= -0.001125 


c = 0.000005 


L = 0.1m 


For carbon steel 

Base temperature is taken as 323K. 


A _ 

2k\a 

ki is the thermal conductivity at base temperature. 
See Fig. 3.13 . 


(3.5) 
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Comparison of T Vs X for constant conductivity and variable conductivity for rectangular fin of a^jper 



in os loooinr^iQcoin 

o C9 o o o 


ajnjBjaclujai iBuoisuaaijpuoN 

Figure 3.12; Comparison of T Vs X for constant and variable thermal con- 
ductivity for copper for parabolic profile 
’ck’ indicates constant conductivity. ’vk’ indicates variable 
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Nondimensional distance 
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Figure 3.13: Comparison of Efficiency Vs v'A for constant conductivity and 
variable conductivity for parabolic fin 

’0.0’ indicates = 0.0, ’0.5 ’indicates ^ = 0.5, ’0.7’ indicates ^ = 0.7 
’0.9 ’indicates = 0.9 , 

For variable conductivity material ’O.Ov’ indicates ^ = 0.0, 

’0.5v’ indicates ^ = 0.5, ’0.7v’ indicates = 0.7, ’0.9v’ indicates ^ = 0.9 . 
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Chapter 4 


OPTIMIZATION OF HEAT 
TRANSFER FROM THE FIN 


4. 1 BASIC OBJECTIVES OF OPTIMIZATION 

Radiating fins are used in spacecraft and the space vehicles for rejecting 
on-board waste heat to deep space . Since weight is at premium , it is 
important to have minimum fin mass in these application or maximization 
of heat transfer per unit weight of the fin. 

4.1.1 DEFINITION OF OBJECTIVE FUNCTION 


The objective function is the heat transfer from the fin per unit mass. 



V is the volume of the fin. 
For rectangular fin 


V = WL (4.2) 

For trapezoidal fin 

(4.3) 

For parabolic fin 

V = a£ + 6y + cy (4.4) 

For all cases , the base temperature (to) is taken as 323K. 

4.1.2 DESIGN VARIABLES 

The design variables are the dimension of the fin for rectangular fin W and L 
,for trapezoidal fin Wi, TF 2 and L . The design variables for parabolic profile 
are a , b , c and L . The heat transfer per unit mass is also function of the 
thermal conductivity . 


4.1.3 LIMITATION AND CONSTRAINTS 

The constraint is the volume of the fin . The limitation is the melting point 
of the fin material. 
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4.2 OPTIMIZATION METHOD 


EXHAUSTIVE SEARCH:-The method consists of calculating the value of 
objective function at the values of variable that are spaced uniformly through 
out the domain of interest. The optimum is obtained by eliminating regions 
not containing the maximum. 

The interval of uncertainty for 1** 2*“* variables are /i, h In- 



Iqi 

Ni + 1 


( 4 . 5 ) 


l2 = 


Io2 

N2 + I 


( 4 . 6 ) 


In = 


Iqn 
Nn + 1 


( 4 . 7 ) 


Where Iqi , I02 Ion are the range of interest for 1** , 2"^ AT*'* 

variables. 

Where number of observation are Ni x N 2 x x Nn ■ 
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4.3 RESULTS FOR ALL THREE TYPE OF 
FINS FOR DIFFERENT MATERIALS 

4.3.1 FOR RECTANGULAR FIN OPTIMUM DIMEN- 
SIONS FOR DIFFERENT MATERIAL 


The calculated optimum dimensions of fin 
RANGE OF INTEREST 
Range of L ( 0.075m, 0.30m) 

Range of IF (0.0001m ,0.0005m) 
Volume of the fin =2.8125 x 10“®7n® 


Type of material 


Length (m) 

Heat transfer(W/kg) 

Aluminium 

0.000163 

0.172546 

396.242 



0.0988567 

78.51465 

Copper 

0.000133 

0.2114662 

147.0347 


Table 4.1: Optimum dimension for Rectangular fins 


4.3.2 FOR TRAPEZOIDAL FIN 

The calculated optimum dimensions of fin 
RANGE OF INTEREST 
Range of L ( 0.075m, 0.30m)) 

Range of Wi (0.0001m, 0.0005m) 

Range of W 2 (0.0001m, 0.0005m) 

Volume of the fin =2.8125 x 10“®m® 
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Type of material 

Wi 

W2 

Length 

Heat transfer(W/kg) 

Aluminium 

0.000267 

0.0001 

0.153270 

556.62384 

Carbon steel 

1 0.000468 

0.0001 

0.0990316 

124.3621 

Copper 

0.000219 

0.0001 

0.176333 

193.5711 


Table 4.2: Optimum dimension for Trapezoidal fins 


4.3.3 FOR PARABOLIC FIN 


The calculated optimum dimensions of fin 
RANGE OF INTEREST 
Range of T ( 0.075m, 0.30m) 

Range of a (0.0001 ,0.0005) 

Range of b (-0.001 , 0.0) 

Range of c (-0.02 ,0.02) 

Volume of the fin =2.8125 x 10“®m^ 


Type of material 

a 

h 

C 

L(m) 

Fin width at tip (m) 

gl(W/kg) 

Aluminium 

0.000121 

-0.001 

0.0029 

0.25 

0.000104 

725.47 

Carbon steel 

0.0002207 

-0.001 

-0.00837 

0.095 

0.0001004 

113.226 

Copper 

0.000184 




0.001 

250.94 

1 


Table 4.3: Optimum dimension for parabolic fins 


4.4 CHOICE OF MATERIAL FOR EACH TYPE 
OF FINS 


Aluminium is the most suitable material for each type of fin because it gives 
maximum heat transfer per unit mass . Which is expected also due to its 
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high conductivity and least density. 
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Chapter 5 


CONCLUSION 


A heat transfer analysis of radiating fins of various materials and various ge- 
ometries is presented. Both constant and variable conductivity fin materials 
have been studied. The results for rectangular fin show excellent agreement 
with an earlier NASA work. The second part of this thesis deals with finding 
dimensions corresponding to the maximum heat transfer per unit weight of 
the fin material having the temperature dependent conductivity. An exhaus- 
tive search method has been used. It is found that aluminium is the most 
suitable material. Also, trapezoidal and parabolic fins turned out to better 
than rectangular fins for efficient heat transfer. 
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Chapter 6 
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Appendix A 


Five point forward differencing 


Five point forward differential for calculating temperature gradient for cal- 
culation of heat transfer. 


, dT . — 25Tj 48Ti+i — 36Ti.f2 “I" 16Tj^3 — 

W''* “ 12AX 
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Appendix B 


Properties of material 

Thermal conductivity k for alumintun = (207.0615 — 0.0054894t + 2.9267 x 
w/m.K 

Thermal conductivity k for carbon steel = (40.882 — 0.0141< ) W/m.K 
Thermal conductivity k for copper =( 390.9975 — 0.0511051t ) W/m.K 
t is in C. 

The above thermal conductivity values are applicable in the range of 
0 - 400 C 

Density of aluminium = 2707 ^ 

Density of carbon steel = 7833^ 

3 

Density of copper = 8954^ 

€ = 0.50 
a = 0.10 

Irradiation G =315.4 ^ 

rnr 
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Appendix C 


Computer Programs 

C.l For Rectangular Fin 

c THE PROGRAM FOR RECTANGULAR FIN 

c THE OUTPUT CAN BE SEEN BY VIEWING ’OUTP’FILE. 
c THIS PROGRAM MUST BE COMPILED ALONGWITH tdma.f. 
c at REPRESENT THE NON DIMENSIONAL TEMPERATURE 
c att REPRESENT THE NON DIMENSIONAL TEMPERATURE IN 
PREVIOUS ITERATION 

DIMENSION a(121),b(121),c(121),d(121),at(121),ak(121),att(121) 

OPEN(UNIT=76,FILE=’OUTP’) 

read *, n 

c*******PROPERTIES AND DIMENSION OF THE FIN ********************* 
e = 0.5 
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alpha = 0.1 

to = 323.14 

al = 0.075 

al = 40.882 

bl = 0.0141 

w = 0.000375 

g = 315.4 

sigma = 5.67e-8 

ts = (alpha*g/e*sigma)**0.25 

const =(w/ (e*sigma*(al**2)*(t0**3))) 

do 10 i= l,n 

c*******CALCULATION OF THERMAL CONDUCTIVITY AT BASE 
TEMPERATURE FOR EACH GRID 
ak(i) = al-(bl*(t0-273.14)) 

c*******ASSUMING THE PREVIOUS ITERATION NON DIMENSIONAL i 
TEMP. FOR EACH GRID 
att(i) = 1.0 

10 continue 

45 ats = ts/tO * 

dx = 1.0/float(n-l) 

^5ie*5ic5|c5|c)fc5iejy|-^jj^ DIAGONAL COFFICIENT OF TDMA ************************4 

b(i)= 1. 

do 30 i= 2,n ; 

i 

I 
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b(i) = ((-2*const*a]£{i)/dx**2)-att(i)**3) 

30 continue 

c*******SUPER-DIAGONAL COFFICIENT OF TDMA *************=(=******** 
c(l) = 0 
do 40 i= 2,n-l 

c(i) =(const*ak(i)/dx**2) + const*(ak(i+l) - ak(i-l))/4.*dx**2 
40 continue 

COFFICIENT OF TDMA *************************= 

a(2) = 0.0 
do 50 i= 3,n-l 

a(i) = ((const*ak(i)/dx**2) - (const*(ak(i+l)-ak(i-l))/2.*dx)) 

50 continue 

a(n) = 2*const*ak(n)/dx**2 

c*******COFFICIENT OF RIGHT HAND COLUMN VECTOR 
d(l) = 1. 

d(2) = -1.0*((const*ak(2)/dx**2)-(const*(ak(3)-ak(2))/4.*dx**2)) 

ats**4 

do 60 i = 3,n 
d(i) = -ats**4 
60 continue 

CALL TDMA(l,n,A,B,C,D,AT) 
do i=l,n 
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c*******CALCULATION OF THERMAL CONDUCTIVITY AT NEW 
TEMP. FOR EACH GRID POINT 

ak(i) = al-(bl*(t0*at(i)-273.14)) 

err=abs(at(i)-att(i)) 

att(i)=at(i) 

if (err.le.O.OOl) go to 55 
go to 45 
55 enddo 

c****** CALCULATION OF TEMPERATURE GRADIENT, HEAT TRANS- 
FER EFFICIENCY ********** 

q = (ak(l)*w*t0*(at(3)-4.*at(2)-t- 3.*at(l))/(2.*dx*al)) 
effeciency = q/((e*sigma*(t0)**4-alpha*g)*al) 

WRITE(76,9) 

9 FORMAT(4x,’THE SOLUTION BY TDMA’/) 

WRITE(76,2) (at (i) ,i= 1 ,n) 

2 FORMAT(2x,5F6.3) 
write(76,8) 

8 forinat(4x,’HEAT LOSS BY THE FIN’/) 
write(76,3) q 

3 format(2x,f9.3) 
write(76,7) 

7 fonnat(4x, ’effeciency of fin’) 

write(76,l) effeciency 
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1 foniiat(2x,f6.3) 
STOP 
END 


c THE NAME OF THE PROGRAM:tdma.f 
c SUBROUTINE FOR SOLVING A SYSTEM OF LINEAR 
c SIMULTANEOUS ALGEBRAIC EQUATIONS HAVING 
c A TRIDIAGONA COEFFICIENT MATRIX, 
c THE EQUATIONS ARE NUMBERED FROM I THROUGH N, 
c AND THEIR SUB-DIAGONAL, DIAGONAL AND SUPER-DIAGONAL 
c ELEMENTS ARE STORED IN THE ARRAYS A,B,AND C.THE 
c RIGHT-HAND COLUMN VECTOR ELEMENTS ARE STORED 
cINTHE ARRAYD.THE COMPUTED SOLUTION VECTORX(I)....X(N) 
c IS STORED IN THE ARRAY X. 

SUBROUTINE TDMA(I,N,A,B,C,D,X) 

DIMENSION A(1),B(1),C(1),D(1),X(1),BETA(101),GAMMA(101) 

^ COMPUTE INTERMEDIATE ARRAYS BETA AND GAMMA 


BETA(I)=B(I) 

GAMMA(I)=D(I)/BETA(I) 

11 = 1+1 
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DO 1 J=I1,N 

BETA{J)=B(J)-A(J)*C(J-1)/BETA(J-1) 

1 GAMMA(J)=(D(J)-A(J)*GAMMA(J-1))/BETA(J) 

c ****** COMPUTE THE SOLUTION VECTOR X ****** 
X(N)=GAMMA(N) 

N1=N-I 
DO 2 K=1,N1 
J=N-K 

2 X(J)=GAMMA(J)-C(J)*X(J+1) /BETA(J) 

END 

C.2 For Trapezoidal Fin 

c THE PROGRAM FOR TRAPEZOIDAL FIN 

c THE OUTPUT CAN BE SEEN BY VIEWING ’OUTP’FILE. 
c THIS PROGRAM MUST BE COMPILED ALONGWITH tdma.f. 
c at(i) REPRESENT THE NON DIMENSIONAL TEMPERATURE 
c ak(i) REPRESENT THE THERMAL CONDUCTIVITY AT ANY 
GRID POINT 

c att(i) REPRESENT THE NON DIMENSIONAL TEMP. IN PREVI- 
OUS ITERATION 

DIMENSION a(121),b(121),c(121),d(121),at(121),ak(121),att(121),x(50) 
DIMENSION cl (100) ,c2(100) ,dl (1 11) ,d2(100) ,al (100) ,a2(100) ,s(100) 
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OPEN(UNIT=76,FILE=’OUTP’) 
read *,n 

c ****** PROPERTIES AND DIMENSION OF THE FiN********=t:=«=f:=c=f!****=c*=i=>f=***= 
e = 0.500 
wl = 0.000375 
alpha = 0.10 
V = 2.8125e-05 
ul = 390.9975 
u2 = 0.0511051 
to = 323. 
dens = 8954 
w2 =0.000100 
al = (2.*v/(wl+w2)) 
g = 315.4 
sigma = 5.67e-8 
ts = (alpha*g/e*sigma)**0.25 
const =(l./(e*sigma*(al**2)*(t0**3))) 
dx = l./float(n-l) 
do 10 i= l,n 
s(l) = wl 

s(i+l) = s(i)-(wl-w2)*dx 

c******CALCULATION OF THERMAL CONDUCTIVITY AT BASE 
TEMPERATURE FOR EACH GRID 
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aJ£(i) = ul-(u2*(t0-273.14)) 

c*******ASSUMING THE PREVIOUS ITERATION NON DIMENSIONAL 
TEMPERATURE FOR EACH GRID 
att(i) = 1.0 
10 continue 
45 ats = ts/tO 

c*******MAIN DIAGONAL COFFICIENT OF TDMA =t:****=t=**********=f=*****=i==f=** 
b(i)= 1. 
do 30 i= 2,n 

b(i) = ((-2*s(i)*const*ak(i)/dx**2)-att(i)**3) 

30 continue 

COFFICIENT OF TDMA ***=i'*******xc=c****=t;*=f:**=(=*=c* 

c(l) = 0.0 

cl(i) =(const*ak(i)*s(i)/dx**2) + s(i)*const*(ak(i+l) - ak(i- 
l))/4.*dx**2 

c2(i) = const*ak(i)*(wl-w2)/2*dx 
c(i) = cl(i)-c2(i) 

40 continue 

^:<c5ic5{c5jc:ic*3Kg|jg_^jQjQQQp^^j^ COFFICIENT OF TDj^^****************************^ 

a(2) = 0.0 
do 50 i= 3,n-l 

al(i) = s(i)*((const*ak(i)/dx**2) - (const*(ak(i+l)-ak(i“l))/2.*dx)) 
a2(i) = const*ak(i)*(wl-w2)/2.*dx 
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a(i) = al(i)+a2(i) 

50 continue 

a(n) = 2*const*ak(n)*s(n)/dx**2 

c*******COFFICIENT OF RIGHT HAND COLUMN VECTOR****************** 
d(l) = 1. 

dl(2) =-1.0*((const*ak(2)/dx**2)-(const*(ak(3)-ak(2))/4.*dx**2))*s(2) 
d2(2) = -(const*ak(2)*(wl-w2)/2 *dx)-ats**4 
d(2) = dl(2)+d2(2) 
do 60 i = 3,n 
d(i) = -ats**4 
60 continue 

CALL TDMA(l,n,A,B,C,D,AT) 
do i=l,n 

c****** CALCULATION OF THERMAL CONDUCTIVITY AT NEW 
TEMPERATURE FOR EACH GRID POINT 
ak(i) = ul-(u2*(t0*at(i)-273.14)) 
err=abs (at (i)-att (i) ) 
att(i)=at(i) 

if (err.le.O.OOOl) go to 55 
go to 45 
55 enddo 

^♦*:(=****q^LqULATION of temperature GRADIENT, HEAT TRANS- 
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FER AND efficiency********* 

anu=-(-25.*at(l)+48.*at(2)-36.*at(3)+16.*at(4)-3*at(5)) 

den = 12 *dx*al 

div = anu/den 

ql = ak(l)*wl*tO*div 

x(l) = 0.0 

do i=l,n-l 

x(i+l)=x(i)+dx 

enddo 

effeciency = ql/((e*sigma*(t0)**4r-alpha*g)*al) 

q = ql/(dens*v) 

print *, wl,w2,al,q, effeciency 

WRITE(76,9) 

9 FORMAT(4x,’THE SOLUTION BY TDMA’/) 

WRITE (76,2) (x(i) ,at (i) ,i=l ,n) 

2 FORMAT(2x,2F6.3) 
write(76,8) 

8 format(4x,’HEAT LOSS BY THE FIN’/) 
write(76,3) q 

3 forniat(2x,£15.3) 
write(76,7) 

7 format(4x, ’effeciency of fin’) 
write(76,l) effeciency 
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1 format (2x,jF9.3) 

STOP 

END 

c THE NAME OF THE PROGRAM:tdma.f 
c SUBROUTINE FOR SOLVING A SYSTEM OF LINEAR 
c SIMULTANEOUS ALGEBRAIC EQUATIONS HAVING 
c A TRIDIAGONA COEFFICIENT MATRIX, 
c THE EQUATIONS ARE NUMBERED FROM I THROUGH N, 
c AND THEIR SUB-DIAGONAL, DIAGONAL AND SUPER-DIAGONAL 
c ELEMENTS ARE STORED IN THE ARRAYS A,B,AND C.THE 
c RIGHT-HAND COLUMN VECTOR ELEMENTS ARE STORED 
c IN THE ARRAY D.THE COMPUTED SOLUTION VECTOR X(I)....X(N) 
c IS STORED IN THE ARRAY X. 

SUBROUTINE TDMA(I,N,A,B,C,D,X) 

DIMENSION A(1),B(1),C(1),D(1),X(1),BETA(101),GAMMA(101) 
c ****** COMPUTE INTERMEDIATE ARRAYS BETA AND GAMMA 

BETA(I)=B(I) 

GAMMA(I)=D(I) /BETA(I) 

11 = 1+1 
DO 1 J=I1,N 

BETA(J)=B(J)-A(J)*C(J-1)/BETA(J-1) 
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1 GAMMA(J)=(D(J)-A(J)*GAMMA(J-1)) /BETA(J) 

C ****** COMPUTE THE SOLUTION VECTOR 

X(N)=GAMMA(N) 

N1=N-I 
DO 2 K=1,N1 
J=N-K 

2 X( J)=GAMMA(J)-C(J)*X( J+1) /BETA( J) 

END 

C.3 For Parabolic Fin 

c THE PROGRAM FOR PARABOLIC FIN 

c THE OUTPUT CAN BE SEEN BY VIEWING ’OUTP’FILE. 
c THIS PROGRAM MUST BE COMPILED ALONGWITH tdmai. 
c at REPRESENT THE NON DIMENSIONAL TEMPERATURE 
c att REPRESENT THE NON DIMENSIONAL TEMPERATURE IN 

PREVIOUS ITERATION 

DIMENSION a(121),b(121),c(121),d(121),at(121),ak(121),att(121) 

DIMENSION cl(100),c2(100),dl(lll),d2(100),al(100),a2(100),s(100) 

DIMENSION x(121),z(121),xx(50) 

OPEN(UNIT=76,FILE=’OUT0 
read *, vO 

PROPERTIES AND DIMENSION OF THE FIN **:»==•=******************> 
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n = 25 
e = 0.500 
alpha = 0.10 
vol = 2.8125e-05 
al = 0.229 
vl = -0.0001000 

v2 = (l./(al*al*al))*((vol/2.)-(v0*al)-(vl*al*al/2.)) 

to = 823. 

ul = 390.9975 

u2 = 0.0511051 

dens = 8954 

wl = 2.*v0 

w2 = 2.0*{v0-(-vl*al-l-v2*al*al) 

g = 315.4 

sigma = 5.67e-8 

ts = (alpha*g/e*sigma)**0.25 

const =(!./ (e*sigma*(al**2)*(t0**3))) 

dx = 1. /float (n-1) 

do 10 i= l,n 

x(l) = 0.0 

x(i-|-l)=x(i)-+-dx 

s(i) = vO + (vl*x(i)*al)-h(v2*((x(i)*al)**2)) 
z(i) = vl*al4-2.*v2*x(i)*al*al 
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c ak(i) = ul 

c*******CALCULATION OF THERMAL CONDUCTIVITY AT BASE 
TEMPERATURE FOR EACH GRID POINT 
ak(i) = ul-(u2*(t0-273.14)) 

c******ASSUMING THE PREVIOUS ITERATION NON DIMENSIONAL 
TEMPERATURE FOR EACH GRID 


att(i) = 1.0 
10 continue 
45 ats = ts/tO 

c*******MAIN DIAGONAL COFFICIENT OFTDMA *************************’ 
b(l)= 1. 


do 30 i= 2,n 

b(i) = ((-4.*s(i)*const*ak(i)/dx**2)-att(i)**3) 


30 continue 


,,*,««« SUPER-DIAGONAL COFFICIENT OF tdIMA***'***************’****'"'’ 


c(l) = 0 


do 40 i= 2,n-l 

cl(i) =(ak(i)*s(i)/dx**2) + s(i)*(ak(i+l) - ak(i-l))/4.*dx**2 
c2(i) = ak(i)*z(i)/2*dx 

c(i) = (cl(i)+c2(i))*2.*const i 

! 

40 continue 

^ ****=f=*SUBDIOGONAL COFFICIENT OF TDMA ***************************** 


a(2) = 0.0 
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do 50 i= 3,n-l 

al(i) = s(i)*((const*ak(i)/dx**2)-(const*(ak(i+l)-alc(i-l))/4 *dx*dx)) 
a2(i) = const*ak(i)*z(i)/2.*dx 
a(i) = (al(i)-a2(i))*2. 

50 continue 

a(n) = 4*const*ak(n)*s(n)/dx**2 

c*******COFFICIENT OF RIGHT HAND COLUMN VECTOR ***=t=****=t=*=f=*=f:***=^ 

d(i) = 1. 

dl(2) =-1.0*((const*ak{2)/dx**2)-(const*(ak(3)-ak(2))/4.*dx**2))*s(2) 
d2(2) = 2.*(const*ak(2)*z(i)/2.*dx)-ats**4 
d(2) - 2.*dl(2)+d2(2) 
do 60 i = 3,n 
d(i) = -ats**4 
60 continue 

CALL TDMA(l,n,A,B,C,D,AT) 
do i=l,n 

c*******CALCULATION OF THERMAL CONDUCTIVITY AT NEW 
TEMP. FOR EACH GRID POINT 

ak(i) = ul - u2*(at(i)*t0-273.14) 
err=abs(at(i)-att (i)) 
att(i)=at(i) 

if (err.le.0.0001) go to 55 
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go to 45 
55 enddo 


c*******CALCULATION OF TEMPERATURE GRADIENT, HEAT TRANS- 
FER AND EFFICIENCY 

aiiu=-(-25*at(l)+48*at(2)-36*at(3)+16.*at(4)-3*at(5)) 

den = 12.*dx*al 

div = anu/den 

ql = ak(l)*wl*tO*div 

eff = ql/((e*sigma*(t0)**4-alpha*g)*al) 

q = ql/(dens*vol) 

xx(l) = 0.0 

do i = 1, n-1 

xx(i-j-l) = xx(i) + dx 

enddo 

print 25, vO ,vl,v2,al,w2,q,eff 
25 fonnat(/7(fll.6)) 

WRITE(76,9) 

9 FORMAT(4x,’THE SOLUTION BY TDMA’/) 

WRITE(76,2)(xx(i), at(i),i=l,n) 

2 FORMAT(2x,2F6.3) 
write(76,8) 

8 format(4x,’HEAT LOSS BY THE FIN’/) 
write(76,3) q 
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3 format(2x,fl5.3) 
write (76, 7) 

7 format (4x,’efFeciency of fin’) 
write(76,l) effeciency 
1 format(2x,f9.3) 

STOP 

END 

c THE NAME OF THE PROGRAM:tdma.f 
c SUBROUTINE FOR SOLVING A SYSTEM OF LINEAR 
c SIMULTANEOUS ALGEBRAIC EQUATIONS HAVING 
c A TRIDIAGONA COEFFICIENT MATRIX, 
c THE EQUATIONS ARE NUMBERED FROM I THROUGH N, 
c AND THEIR SUB-DIAGONAL, DIAGONAL AND SUPER-DIAGONAL 
c ELEMENTS ARE STORED IN THE ARRAYS A,B,AND C.THE 
c RIGHT-HAND COLUMN VECTOR ELEMENTS ARE STORED 
c IN THE ARRAY D.THE COMPUTED SOLUTION VECTOR X(I)....X(N) 
c IS STORED IN THE ARRAY X. 

SUBROUTINE TDMA(I,N,A,B,C,D,X) 

DIMENSION A(1),B(1),C(1),D(1),X(1),BETA(101),GAMMA(101) 
c ****** COMPUTE INTERMEDIATE ARRAYS BETA AND GAMMA 

BETA(I)=B(I) 
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GAMMA(I)=D(I) /BETA(I) 

11 = 1+1 
DO 1 J=I1,N 

BETA(J)=B(J)-A(J)*C(J-1)/BETA(J-1) 

1 GAMMA(J)=(D(J)-A(J)*GAMMA(J-1))/BETA(J) 

c ****** COMPUTE THE SOLUTION VECTOR X ****** 
X(N)=GAMMA(N) 

N1=N-I 
DO 2 K=1,N1 
J=N-K 

2 X(J)=GAMMA(J)-C(J)*X{J+1) /BETA(J) 

END 
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